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Abstract 

System-size dependence of dihadron azimuthal correlations in ultra-relativistic heavy ion colli- 
sion is simulated by a multi-phase transport model. The structure of correlation functions and 
yields of associated particles show clear participant path -length dependences in collision systems 
with a partonic phase. The splitting parameter and root-mean-square width of away-side correla- 
tion functions increase with collision system size from 14 N+ 14 N to 197 Au+ 197 Au collisions. The 
double-peak structure of away-side correlation functions can only be formed in sufficient "large" 
collision systems under partonic phase. The contrast between the results with partonic phase and 
with hadron gas could suggest some hints to study onset of deconfinment. 
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1. Introduction 

Quantum Chromodynamics (QCD) calculation predicted an exotic quark-gluon matter in 
QCD phase diagram yj] may be created in the early stage of heavy ion collisions at ultra- 
relativistic energy fl. For mapping the QCD phase diagram and locating QCD phase boundary 
and critical point yfl, one needs to find a way to vary temperature (T) as well as chemical po- 
tential (/is). The NA61 collaboration and NA49-future collaboration |4j] suggested that it can be 
achieved via a systematic energy (E) and system-size (A) (E - A) scan. 

Jet quenching phenomenon has been theoretically predicted and experimentally ob- 
served [6J. So far, dihadron azimuthal correlations have been demonstrated as a good method to 
reconstruct particle and energy distribution induced by the quenched jet. In experiment, a double- 
peak structure was found on the away side of dihadron azimuthal correlation functions El 
and the indication of conical emission of charged hadrons was reported by the STAR collab- 
oration [10]. The centrality and transverse momentum dependences of double-peak structure 
of away-side correlation functions were experimentally investigated by RHIC-BNL HEl and 
theoretically simulated in Ref . 11211 . 

These interesting phenomena have attracted some theorists to explain the physical mecha- 
nisms for the origin of the double-peak structure. These mechanisms include a Cherenkov-like 
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gluon radiation model 1 131. medium-induced gluon bremsstrahlung radiation U3. Q, shock 
wave model in hydrodynamic equations [16], waking the colored plasma and sonic Mach cones JT 
sonic booms and diffusion wakes in therma l gauge -string duality 1 181. jet deflection lfl9tl and 
strong parton cascade mechanism and so on (12, 2(1 21 , 22, 23]. Recently, Gyulassy and his col- 
laborators suggest that the conical emission can stem from universal flow-driven mechanism |24 ] . 
In Ref. l25ll . the shock wave phenomena are discussed in viscous fluid dynamics and kinetic the- 
ory. Renk and Neufeld II26I1 presented their systematic study of dependence of the Mach cone 
signal on the energy deposition into the medium in linearized hydrodynamics. The path-length 
effect on the correlations relative to the reaction plane are studied, respectively, by Jia et al. in 
a simple model [27] and by Ma et al. in AMPT model f28ll . In the shock wave model 1 16], the 
emission angle relative to jet is calculated to be about 1.23 rad for QGR 1.11 rad for hadronic 
gas and zero for mixed phase. The gluon radiation mechanism for the double-peak structure IU3I1 
suggests that the more energetic the jet, the smaller the emission angle. In Ref. i29tl . it is sug- 
gested that the suppression or even disappearance of Mach cone at the QCD critical point due 
to the attenuation of the sound mode. While it remains unclear what the main mechanism for 
the emergence of the double peak structure is, all the experimental and theoretical works sugges t 
that it should depend on the nature of the hot and dense matter created in the collisions DOiblll . 
In this paper, we study the properties of hot-dense matter produced by different system size by 
investigating the system-size dependence of dihadron azimuthal correlations. 

In this paper, we present participant path-length, defined as v = 2Nbi„/N par t l32ll (Nbin and 
Npan are the number of binary collision and participants, respectively), dependence of the double- 
peak structure of away-side correlation function in the most central collisions (0-10%). The 
structure of away-side correlation function changes near 40 Ca + 40 Ca collisions at yjs^N = 200 
GeV in central collisions (0-10%). The results show obvious degree of freedom dependence in 
the system with a partonic phase or with a pure hadron gas 113 311 . which implies information of 
the onset of deconfinement. 



2. Model and analysis method 



In this work, a multi -phase transport model (AMPT) [34], which is a hybrid dynamic 
model, is employed to study dihadron azimuthal correlations. It includes four main compo- 
nents to describe the physics in relativistic heavy ion collisions: 1) the initial conditions from 
HIJING model f35ll . 2) partonic interactions modeled by a Parton Cascade model (ZPC) lf36Tl . 
3) hadronization (discussed later), 4) hadronic rescattering simulated by A Relativistic Transport 
(ART) model 13711 . Excited strings from HIJING are melted into partons in the AMPT version 
with string melting mechanism ]38|] (abbr. 'the Melt AMPT version ') and a simple quark coa- 
lescence model is used to combine the partons into hadrons. In the default version of AMPT 
model {39I (abbr. 'the Default AMPT version'), minijet partons are recombined with their parent 
strings when they stop interactions and the resulting strings are converted to hadrons via the Lund 
string fragmentation model ll40ll . The Melt AMPT version undergoes a partonic phase, while a 
pure hadron gas is in the Default AMPT version. Details of the AMPT model can be found in a 
review paper [34] and previous works 1341 38l 41 1. 

The analysis method for dihadron azimuthal correlations is similar to that used in previ- 
ous experiments iBTl lill. which describes the azimuthal correlation between a high pj particle 
(trigger particle) and low pj particles (associated particles). The raw signal can be obtained by 
accumulating pairs of trigger and associated particles into A(f> = <p assoc - ^ t rig distributions in the 
same event. The background which is expected mainly from elliptic flow is simulated by mixing 

2 



event method OIL 1811 . To reconstruct the background, we accumulate pairs of one fake trigger 
particle (high pj) in one event and another fake associated particle (low pj) in another event 
to obtain the A0 distribution as the corresponding background, the centralities of the above two 
events are requested very closed. Then the background is subtracted from raw signal by using "A 
Zero Yield At Minimum" (ZYAM) assumption as that used in experimental analysis 18J] (see our 
detailed analysis in Ref. Il2in . Recently, Wang et al. 14211 discussed the background in the cor- 
relations and presented an analytical form for flow background to jet-like azimuthal correlations 
in a cluster approach. And it is suggested that the collision geometry fluctuations and triangular 
flow should be taken into account in the correlation analysis B43l l44l 14511 . But those go beyond 
our discussion in this paper. 



3. Results and discussions 

3.1. Structure of dihadron correlation function 

The participant path-length, defined as v = 2Nbml Npart 1I32I1 . can describe degree of multiple 



collisions between participants in the early stage of heavy ion collisions and characterize the 
size of the reaction zone. The n'™"°" represents average collision number of partons in the Melt 
AMPT version. The values of v and n pa *, on significantly increase with varying collision system 
(CSYS) from "small" size to "large" size at -\/snn = 200 GeV in the most central collisions (0- 
10%) as shown in Table [TJ From this table, we can see the multiple collisions are more frequent 
in "large" size collision system than in "small" size one. The values of N part , Nbin and v are 



comparable to those from the Glauber Model |]46|] for both 64 Cu + 64 Cu and 197 Au + 197 Au 
collisions. 

When the mixed background which mainly stems from elliptic flow is subtracted from the 
raw dihadron correlation signal taken in the same events, we can get the correlation function. 
Fig.[T]shows dihadron azimuthal correlation functions of different collision systems in the most 
central (0-10%) collisions at yJs^N = 200 GeV. The correlation functions are calculated in the 
kinetic windows, 1 < p™ soc < 3 GeV/c as well as 2.5 < p" 8 < 6 GeV/c and \rj\ < 1. It shows 
that the structure of away-side correlation function changes from the Gaussian-like distribution 
to double-peak structure near 40 Ca + 40 Ca collisions with varying collision system from 14 N + 
14 N to 197 Au + 197 Au collisions in the Melt AMPT version. In this figure, the amplitude of the 
correlation function becomes higher with the increasing of collision system size. The associated 
particles in the Melt AMPT version are more abundant than those in the Default AMPT version. 



Table 1: N parl (CSYS), Nbin(CSYS), v(CSYS) = n'""\csyS) ' n col'°" ' n different collision system at t/snn = 200 GeV 
for centrality 0-10 %, the value in blanket is taken from Glauber Model |46ll . 



CSYS 


14 N + !4 N 


I6q + 16q 


23 Na + 23 Na 


27 Al + 27 Al 


40 Ca + 40 Ca 


64 Cu + 64 Cu 


197 Au + 197 Au 


N, mr ,(CSYS) 


20.78 


24.25 


35.92 


43.61 


65.97 


107.04(99.0) 


343.32 (325.9) 


N hi „(CSYS) 


19.63 


23.69 


41.01 


54.34 


91.15 


179.98 (188.8) 


914.71 (939.4) 


v(CSYS) 


1.89 


1.95 


2.28 


2.49 


2.76 


3.36 (3.8) 


5.33 (5.7) 


'C°" 


1.31 


1.44 


1.93 


2.23 


2.79 


3.80 


7.24 



The Default AMPT version is used to compare the properties of the double-peak structure 
in partonic phase and in hadron gas. For investigating the properties of collision system-size 
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Figure 1: Dihadron azimuthal coiTelation functions for different collision systems in centrality 0-10% at ^/Fnn = 200 
GeV; Kinetic windows: I < p a T sso < 3 GeV/c, 2.5 < p"' g < 6 GeV/c, M <1; the shadowing area from the STAR data |7j]. 

dependences of away-side dihadron azimuthal correlations, we extract the associated particle 
yield N™?°f , splitting parameter (D) (half distance between double peaks on the away side) and 
Root Mean Square Width (A<f> rms ) of away-side associated particles, which will be discussed in 
the following sections, respectively. 



3.2. Yield of associated particles 
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Figure 2: Yield of associated particles on away-side correlation functions, 
Melt/Default AMPT version in the centrality of 0-10% at ^fs^ = 200 GeV. 



N"' 



f, as a function of v(CSYS) by the 



The v(CS YS) dependence of N™™° is shown in Fig.|2]from the Melt/ Default AMPT version, 
respectively. It presents a significant increasing trend of N" t ^ with varying the collision system 



from 14 N + 14 N to 197 Au + 197 Au collisions in the most central collisions (0-10%) at 

4 



105 200 GeV in the Melt AMPT version. The Default AMPT version, with a hadronic gas, does not 

ioe result in a rapid increasing dependence trend. In the Melt AMPT version, the dependence trend 

107 indicates the jet correlation information can be inherited by and transmitted to more particles in 

108 a partonic phase than in a hadronic gas, especially in "large" size collision system. Furthermore 
loo it implies that the interaction strength in "large" size collision system is more significant than 
no that in "small" size collision system, and while strong parton cascade plays a dominant role 
in to push more away-side associated particles in the Melt AMPT version. It is interesting that 

112 the increasing slope of N° S J°^ vs v{CSYS) from the linear fitting in the Melt AMPT version is 

113 quicker above 40 Ca + 40 Ca collision system, where clear double-peak structure emerges, than 

114 that in small systems. This property indicates the double-peak (Mach-like) structure can enhance 

115 associated particles yields of jet correlations partially. 

no 3.3. A(p rms and splitting parameter on the away side 

117 Root Mean Square Width (A<p rms ) of away-side correlation function is defined as 

2 (t4-M> m )H\IN,rig){dNldt4) 

away 

E (l/N lrig )(dN/dAcf>) ' 

away 

us where A<f> m is the mean A<p of away-side correlation function and it approximates to n. A(p rms 

no can describe the diffusion degree of the associated particles relative to the direction of back jet. 

120 The v(CSYS) dependences of A<p rms in the Melt/ Default AMPT version are shown in Fig. [3] 

121 respectively. A0„„ s from the Melt AMPT version are consistent with PHENIX data for 

122 Cu + Cu and Au + Au collisions. A(f> rms increases from 14 N + 14 N collisions to 197 Au + 197 Au 

123 collisions in the Melt AMPT version as well as the Default AMPT version, but the increasing 

124 trend is not so quick in the later, especially for systems larger than 40 Ca + 40 Ca. The increas- 
es ing trend of A<f> rms shows broadening of away-side correlation functions with increasing size of 

126 collision system. It indicates that the jet correlation information can reach faraway relative to 

127 direction of jet with changing the collision system from "small" size one to "large" size in a par- 
es tonic phase. It is remarkable that the increasing trend of A<p rms from the linear fitting in the Melt 

129 AMPT version shows two different slope after and before 40 Ca + 40 Ca collision system, where 

130 clear double-peak structure emerges. The double peak (Mach-like) phenomenon from quenched 

131 jet can enhance the diffusion degree of the associated particles relative to back jet in the "large" 

132 size collision system. This suggests that the back jet modification in the medium created in heavy 

133 ion collisions with a partonic phase is more distinct in the "large" size collision system than in 

134 "small" size one. 

135 The splitting parameter (D) is another useful observable to characterize the structure of the 

136 double-peak of away-side correlation function, and further discloses essential of jet modification. 

137 The v(CSYS) dependence of splitting parameters (D) in the Melt/Default AMPT version are 

138 shown in Fig. [4] for the most central collisions at -\/snn = 200 GeV. For lighter systems from 

139 14 N + 14 N to 27 Al + 27 A1 collisions, the splitting parameter is not extracted since there is no 

140 observable double-peak structure of away-side correlation functions. In both simulation cases, 

141 the splitting parameter (D) increases from "small" size collision system to "large" size one, which 

142 indicates that there exists stronger jet-medium interaction in "large" system. It is also remarkable 

143 that the splitting parameter (D) is larger in the Melt AMPT version than that in the Default AMPT 

144 version. The Melt AMPT results are comparable to PHENIX data 19fl_for Cu + Cu and Au + Au 

145 collisions due to effect of parton cascade in the Melt AMPT version 11211 . The parton interaction 
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Figure 3: RMS-Width (A0 r „„) of the away-side correlation functions as a function of v(CSYS ) in the Melt/Defualt 
AMPT version for the centrality 0-10% at ^fsJ^ = 200 GeV; Square from PHENIX data 
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cross section is taken to be 10 mb in this work, which is also reasonable fo r rep roducing elliptic 
flow as well as dihadron azimuthal correlations in the Melt AMPT version l34ll38ll4Ull2ll . This 
calculation reflects that it is necessary to pass a strong partonic stage to reproduce large enough 
double-peak structure as the experimental data demonstrate. 

From these results, it can be concluded that a considerable "large" collision system is neces- 
sary and the strong parton cascade is essential for the formation of the double-peak structure of 
away-side correlation function. An onset of the observable double-peak structure occurs in the 
mass range of 40 Ca + 40 Ca collision. This phenomenon indicates that the correlation is sensitive 
to v(CS YS) and n pa , °", i.e. the correlation depends on the collision system size and the violent 
degree of the partonic interaction in a partonic phase. Different results obtained in a partonic 
phase and a pure hadron gas imply the double-peak structure and jet modification are sensitive to 
the effective degree of freedom of the dense medium created in relativistic heavy ion collisions, 
which can give us some hints of the onset of deconfinement in the system-size viewpoint. 



159 4. Summary 

160 In summary, the present work discusses the collision system-size dependence of dihadron 

161 azimuthal correlations at ysjvw = 200 GeV by a multi-phase transport model. The yields of 

162 associated particles, width of away-side correlation functions and splitting parameter show sig- 

163 nificant system-size dependence. The away-side correlation function becomes more broadening 

164 with the increasing of collision system size and displays the onset of double-peak structure near 

165 40 Ca + 40 Ca collisions after the system undergoes a strong partonic transport stage. These results 

166 also present the degree of freedom dependence, which might be related to onset of deconfine- 

167 ment. We would remark that these observations do not assume any dynamical mechanism for the 
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Figure 4: Splitting parameter (D) of the away-side correlation functions as a function of v(CS YS) in the Melt/Defualt 
AMPT version for the centrality 0-10% at ^fTM = 200 GeV; Square from PHENIX data £J. 



168 formation of the double hump structure. The AMPT model would include collective Mach-like 

169 effects associated to the particles, but it would also include eccentricity fluctuating effects and 

170 triangle flow components. The result is, in this sense, robust. 
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